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We consider a neutrino mass model where all leptonic mixing is induced by a heavy Majorana 
sector through the seesaw type I mechanism, while the Dirac mass matrices are diagonal. Such a 
pattern occurs naturally in grand unified theories. Constraints on the parameters of the models are 
considered and it is shown that a normal neutrino mass hierarchy is preferred. The lightest neutrino 
mass is typically small, leading to nonobservable rates for neutrinoless double beta decay in the 
normal hierarchy case. 
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I. INTRODUCTION 



In the absence of light right-handed neutrinos, neutrino masses can be generated through the seesaw mechanism. 
The mass matrix of the neutrinos can then be written as 

M^ = m T M - 1 mu, (1) 

where tod refers to the matrix of Dirac masses and M is the matrix of Majorana masses. 

The mixing of the neutrinos can be described by the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix Upmns- 
At the 3cr level it can be experimentally determinedpQ to be, 

/ [0.77, 0.86] [0.50,0.63] [0.00,0.22]\ 
V^pmns- [0.22,0.56] [0.44,0.73] [0.57,0.80] . (2) 
\ [0.21, 0.55] [0.40,0.71] [0.59,0.83]/ 

Within these bounds, the PMNS matrix is compatible with the tribimaximal (TBM) pattern[3]: 

Vpmns ~ ^tbm - 71 73 75 • (3) 
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II. OUTLINE OF THE MODEL 



In this paper we assume that all leptonic mixing originates from a heavy Majorana sector while the Dirac mass 
matrices of the neutrinos and charged leptons are diagonal. 

A model of this type is quite a natural consequence of 5*0(10) grand unified theories (GUTs). With all quarks 
and leptons unified in a 16 multiplet of S'O(IO), the GUT will cause the mass matrices of the quarks and leptons 
to be very similar. This means that they could at least approximately all be brought into a diagonal form. The 
experimentally obvious differences between the CKM quark mixing matrix and the PMNS lepton mixing matrix have 
to be explained by an additional mechanism, which can be the seesaw mechanism used in this paper. 

With respect to the seesaw formula Eq. [I] the model investigated here can be written as 

/m? \ 

m D = ro? , (4) 
V m°/ 
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and M € R 3x3 symmetric and arbitrary. After applying Eq. [T] the mass matrix in the flavor basis is given by 
f hl 22 M 33 - (M 23 ) 2 j (m?) 2 (M 13 M 23 - M 33 M 12 )m?m$ (M 12 M 23 - M 22 M 13 )mfm^\ 



(M 13 M 23 -M 33 M 12 )mYm$ M n M 33 - (M 13 ) 2 (m%) 2 (M 12 M 13 - Af 11 M 23 )m 2 3 m 3 3 



v (Mi 2 M 23 - M 22 M 13 )m¥m$ (MuMia - M 1 iM 23 )m 2 3 m 3 3 



M n A/ 22 - (Af 12 ) 2 J (m») j 



(5) 



where the common factor of mass dimension three is given by A 3 = — M 33 (Afi 2 ) 2 + 2Afi 2 Afi 3 Af 23 — Af 22 (Mi 3 ) 2 — 
Mil (Af 23 ) + AfnAf 22 Af 33 . It is obvious that the structure of the matrix depends crucially on the differences of the 
Major ana masses Afy. 

The (1, 1) element of Eq. [5]is the effective mass mpp observed in neutrinoless double beta decays. In order to make 
statements about the neutrino mixing angles and squared mass differences, the mass matrix has to be diagonalized 
using an eigenvalue decomposition, yielding the mixing matrix U. We use the ordering scheme of Ref. [3] in which the 
labels mi and m 2 are assigned to the pair of eigenvalues whose absolute squared mass difference is minimal. Out of 
these two the eigenvalue whose corresponding eigenvector has the smaller modulus in the first component is labeled 
m 2 . The hierarchy of the neutrino masses is then given by the sign of the squared mass differences Am 3i or Am 32 . 
The mixing angles can then be determined using [4] 

arcsin(|i7i 3 |) , 



713 



arctan 



arctan 



'23 



(\Ul2\\ 

\\Uu\J 

\\U33\) 



if U n £ 
else 

if U 33 ? 
else 



(6) 
(7) 

(8) 



If the mass matrix in Eq. [5] is supposed to represent neutrino data, it needs to be able to generate neutrino mixing 
that is close to being tribimaximal[2]. A general class of flavor space mass matrices that leads to tribimaximal mixing 
is given by the following pattern[S]: 



M, 



TBM.flv 



x y y 
y x + v y — v 
y y-v x + v, 



(9) 



where x,y and z are real numbers. It is useful to compare the entries of this matrix to Eq. [5] for the two important 
mass hierarchies: 

1. inverted hierarchy: In this case we approximate an inverted mass hierarchy by two neutrino masses at a 
higher scale in and one neutrino mass set to zero — i.e., the diagonal mass matrix becomes diag(m, in, 0). Eq. [9] 
can then be written as 
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(10) 



which is equivalent to setting x = 1, y = Q and v = — x ji in Eq. [9j Comparing Eq. 10 with the mass matrix of 
Eq. [5] leads to a set of equations whose solutions determine if an inverted neutrino mass hierarchy is possible in 
this model: 
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Trying to solve this set of equations immediately leads to a condition 

m = 0, (12) 

which means that in this approximation it is not possible to generate an inverted neutrino mass hierarchy. 
Translated to a realistic scenario where the facts that tribimaximal mixing is only an approximation and that 
the smaller squared mass difference is not zero are taken into account, one can conclude that in this model the 
inverted mass hierarchy should be strongly suppressed. 

2. normal hierarchy: A normal neutrino mass hierarchy is approximated by two vanishing neutrino masses and 
one neutrino mass at a higher scale rh — i.e., a diagonal mass matrix of diag(0, 0, rh). This leads to a flavor space 
mass matrix of the form 




(13) 

Again, comparing this matrix with Eq. [5]gives a set of equations. Note that in this case the (1, 1) element of the 
matrix is zero instead of rh. This eliminates the suppressive condition of Eq. 12 The rest of the set of equations 
is solvable and just restricts the parameter space of the mass matrix. 



III. NUMERICAL ANALYSIS 



In this section we present the results of a numerical analysis to determine the general allowed ranges for the following 
observables: 

(i) the neutrinoless double beta decay parameter to^/j, given by the (1, l)-entry of Eq. [5j 

(ii) the lightest neutrino mass mo, 

(iii) the neutrino mixing angle $13, as it is the mixing angle with the largest experimental uncertainty. 

The Dirac masses are mostly responsible for the mass eigenvalues of the neutrinos, while the mixing angles are 
dominantly determined by the Majorana mass matrix entries My. 

Because of this, the numerical analysis of each point in the parameter space is performed in two steps: First, random 
starting points of the electroweak scale are chosen for the Dirac masses mj 3 , which are then varied. The Majorana 
parameters My are also chosen randomly at a scale of up to 100 x 10 14 GeV, but are not varied in this step. We 
have ignored the possibility of CP violation here and assigned real valued numbers to all parameters. A x\i function 
for the squared mass differences can be calculated using the diagonalized mass matrix and comparison values from 
a global fit of all experimental data[6j (the two cited values for Am^ refer to the two possible mass hierarchies; the 
uncertainties refer to la and 3<r respectively): 

(Am^) 6XP = 7.59 ± 0.20 (+°^) x 10^ eV 2 , (14) 

CXP= |-2.36 ± 0.11( ± 0.37)xl0- c v 2 
v 31J |+2.46 ±0.12 (±0.37) x 10~ 3 eV 2 v ; 

This Xm function is minimized using a multidimensional minimization algorithm S7 . At this point, if the minimum 
is above the threshold value for Xm, the data point is discarded. For the accepted points, the second step consists of 
calculating the mixing angles [3] . The \\ function for the angles is then analyzed and compared to the data obtained 
from the global fit of experimental results [6]: 



&g = 34.4 ± 1.0 ( ^9) , ^ P - 42.8 + _ A 2 l (T 7 3 7 ) , <ff = 5-6 t 3 2 ° 7 (< 12-5)° . (16) 

Note that the cited value for #13 is based on a global fit including the other angles and mass observables as no direct 
measurement of the angle exists as of now. 

If the x\ value is below the threshold value (for the number of degrees of freedom considered here, it is \\ < 16.8 
for a significance of 0.99), the data point as well as the values of the observables are recorded. 
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The scan covering roughly 16 x 10 6 data points finds 177256 acceptable samples that lead to a normal neutrino mass 
hierarchy and 3596 samples with an inverted neutrino mass hierarchy. This meets the expectation of the inverted 
mass hierarchy being suppressed in the model. The best fit point lies in the regime of normal mass ordering with 
\\ = 0.005. The parameters leading to this best fit are 

Mil = 44 x 10 13 GeV M 22 = 24 x 10 13 GcV M 33 = 89 x 10 13 GeV 

M12 = 97 x 10 13 GeV M 13 = 97 x 10 13 GcV M 23 = 21 x 10 13 GeV (17) 

m? = 46 GcV m° = 140 GeV m° = 130 GeV 

For all acceptable points, the values for the lightest neutrino mass mo are displayed in Fig. [I] for a normal mass 
hierarchy and in Fig. [2] for an inverted neutrino mass hierarchy. 

For both hierarchies, the best fit lightest neutrino mass mo is below 0.005 eV, although larger values of up to 0.02 cV 
are still acceptable for the inverted hierarchy case (Fig. [5]) . Both values lead to summed neutrino masses M„ = m i 
below the current bound of roughly M„ w 0.5 eV[5HT7]. 

As the contribution from mo is negligible and #13 is small, the neutrinoless double beta decay observables for the 
cases of normal or inverted hierarchies are given by (see e.g. Refs. [18-21 and the references therein) 

{y Am] 2 sin 2 (#12) for normal hierarchy 

yj Am% 3 resp. \J Am^ cos (26*12) for inverted hierarchy ' 

where Ami2 <C Am 2 3 and CP conservation has been assumed. The two values given for the case of an inverted 
mass hierarchy stem from the sign ambiguity of the ± sin 2 (#12) term in the sum of masses. Using these formulas, the 
squared mass differences and mixing angles given in the global fit of Ref. [6] lead to m^ rmal ps 0.003 eV for the case 

of a normal mass hierarchy. The two possible values for the inverted mass hierarchy case are rrtaZ ~ 0.02 eV and 

rripZ^ h « 0.05 eV. As the model discussed in this article reproduces the mixing angles and mass parameters of the 
global fit, the predictions for the mass parameter of the neutrinoless double beta decay mpp in that case also hold. 

All angles can be fitted to the experimental data|£y. As there is no direct measurement for 6*13, it is instructive to 
also take a look at the case where #13 is not part of the fit. In that case, only the squared mass differences and the 
angles 6*12 and 6*23 are fitted. The allowed ranges for the angle #13 can then be seen as a prediction of the model. In 
the case of a normal mass hierarchy, there is a clear preference for 6*13 « 0, although values up to and beyond the 
bound of 6^3 P < 0.15(22] are certainly possible. In the case of an inverted mass hierarchy, the preferred value of #13 
is clearly close to a maximal mixing angle, although there are results that satisfy the experimental bound. 

Consequently for the preferred case of normal mass hierarchies the model typically predicts small values for 6*13 
outside of the sensitivity range of Double Chooz [IS EH] , Daya Bay[251 HU and Reno[271 [2S] (see also Refs. 150]). 
For the less favored inverted hierarchy, however, upcoming reactor experiments have good prospects for measuring 

Recently, the T2K experiment [3T] has observed indications of — > v e appearance [32j. The observed number of 
events deviates from the expectation compatible with sin 2 2#i3 = with a signifiance of 2.5er. According to Ref. [32] 
this leads to an allowed range of 0.03(0.04) < sin 2 26>i 3 < 0.28(0.34) for a normal (inverted) mass hierarchy. Even 
though this model prefers a large value for #13 in the inverted hierarchy case, the predictions remain compatible with 
the T2K bounds within the allowed x 2 range for both normal and inverted mass hierarchies. 



IV. CONCLUSIONS 



In this paper, a generic model based on a seesaw type 1 mechanism with diagonal Dirac mass matrices for both the 
charged leptons and the neutrinos has been considered. All contributions to the observed neutrino mixings originate 
from the heavy Majorana masses through a generic Majorana mass matrix that allows for off-diagonal components. 
It has been shown that models of that kind — which are well-motivated by GUTs — generate a small mixing angle 6*13 
naturally and that a normal neutrino mass hierarchy is preferred. 
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FIG. 1. The lightest neutrino mass mo for all models with \a < 16.8 for the case of a normal neutrino mass hierarchy. The 
shade is proportional to the number of hits that lie in the shaded bin. The plot has been divided into 1668 bins on the x axis 
and 35 bins on the y axis. Outside of the boundary line the hit density is less than 10 hits per bin. 
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FIG. 2. The lightest neutrino mass mo for all models with Xa < 16.8 for the case of an inverted neutrino mass hierarchy. The 
shade is proportional to the number of models that lie in the shaded region. The plot has been divided into 689 bins for x axis 
and 35 bins for the y axis. 



